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Abstract. A local turbulence model is developed to study energy cascades in the heliosheath and outer heliosphere (OH) 
based on self-consistent two-dimensional fluid simulations. The model describes a partially ionized magnetofluid OH that 
couples a neutral hydrogen fluid with a plasma primarily through charge-exchange interactions. Charge-exchange interactions 
are ubiquitous in warm heliospheric plasma, and the strength of the interaction depends largely on the relative speed between 
the plasma and the neutral fluid. Unlike small-length scale linear collisional dissipation in a single fluid, charge-exchange 
processes introduce channels that can be effective on a variety of length scales that depend on the neutral and plasma densities, 
temperature, relative velocities, charge-exchange cross section, and the characteristic length scales. We find, from scaling 
arguments and nonlinear coupled fluid simulations, that charge-exchange interactions modify spectral transfer associated with 
large-scale energy-containing eddies. Consequently, the turbulent cascade rate prolongs spectral transfer among inertial range 
turbulent modes. Turbulent spectra associated with the neutral and plasma fluids are therefore steeper than those predicted by 
Kolmogorov's phenomenology. Our work is important in the context of the global heliospheric interaction, the energization 
and transport of cosmic rays, gamma-ray bursts, interstellar density spectra, etc. Furthermore, the plasma-neutral coupling is 
crucial in understanding the energy dissipation mechanism in molecular clouds and star formation processes. 
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1. INTRODUCTION 

In many circumstances, astrophysical and helioshperic 
plasmas are characterized by partially ionized gases in 
which magnetized plasma of protons/electrons and neu- 
tral particles co-exist. Number densities of these compo- 
nents and interaction processes can vary. The plasma and 
the neutral particles especially in the heliospheric plasma 
interact mutually through charge exchange, i.e., a pro- 
ton sufficiently near a neutral hydrogen atom can capture 
its electron, thereby creating a new neutral and a proton. 
This process conserves plasma and neutral densities, but 
not momentum and energy. The charge exchange mean 
free path in the local heliospheric plasma is typically 
about the order of 50 AU [1]. Length-scales smaller than 
the charge exchange mean free path can exhibit turbulent 
motion. The physics of small-scale turbulent motions is 
complex and is important to e.g. heating, the transport of 
cosmic rays, Fermi acceleration, the density spectrum. 
However no self-consistent simulation models attempt 
to directly address the multi-component and multiple- 
scale character of the outer heliosphere and its associated 
turbulence. Here we present a self-consistent plasma- 
neutral turbulence simulation model based on Fourier 
spectral techniques that are widely employed in the stud- 
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ies of neutral and plasma fluids. We concentrate on some 
of the most fundamental aspects of partially ionized he- 
liospheric plasma turbulence. In section 2, we discuss the 
equations of the coupled plasma-neutral fluid. Section 3 
describes results of nonlinear self-consistent fluid simu- 
lations that evolve plasma and neutral fluids in the pres- 
ence of charge exchange forces. The charge exchange 
spectra are discussed in section 4. Conclusions are pre- 
sented in section 5. 



2. PLASMA NEUTRAL COUPLED FLUID 
MODEL 

Our model simulates the partially ionized plasma in a 
local region. The fluid model describing nonlinear turbu- 
lent processes, in the presence of charge exchange, can 
be cast into plasma density (p p ), velocity (Up), magnetic 
field (B), pressure (P p ) components according to the con- 
servative form 



dFp 
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The above set of plasma equations is supplemented by 
V • B = and is coupled self-consistently to the neutral 
density (p„), velocity (V„) and pressure (P n ) through a 
set of hydrodynamic fluid equations, 
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(2) 



where, 





Pn 




p«V„ 


F„ = 


PnVn 


,Q„ = 


Pn * n Vn ~r" y_ j 




<?« 









=s„. p = QmOv^u^p^p,,,?;,:?},) 
_ g£(v„,Up,pp,p„,r„,r p ) 

e " = 2 p " y " + r-T- 

Equations (Q~|i to (O form an entirely self-consistent 
description of the coupled plasma-neutral turbulent 
fluid. The charge-exchange momentum sources in the 
plasma and the neutral fluids, i.e. Eqs. ([TJ and ©, are de- 
scribed respectively by terms Qm(U p ,\ r „,Pp,p n ,T n ,T p ) 
and Qm(V„, U p , p p ,p», T n ,T p ). A swapping of the 
plasma and the neutral fluid velocities in this rep- 
resentation corresponds, for instance, to momentum 
changes (i.e. gain or loss) in the plasma fluid as a 
result of charge exchange with the ISM neutral atoms 
(i.e. QM(V p ,\ n ,p P) p n ,T n ,Tp) in Eq. ££})■ Similarly, 
momentum change in the neutral fluid by virtue of 
charge exchange with the plasma ions is indicated by 
QM(yn,Up,Pp,Pn,T n ,Tp) in Eq. ©. For a complete 
description of charge exchange forces, the readers can 
refer to our work in Ref . 12[] . 



3. Simulations 

A two-dimensional (2D) nonlinear fluid code was de- 
veloped to numerically integrate Eqs. (JT|) to (J2]i. The 
spatial discretization in our code uses a discrete Fourier 
representation of turbulent fluctuations based on a pseu- 
dospectral method |Ht], while we use a Runge Kutta 4 
method for the temporal integration. All the fluctua- 
tions are initialized isotropically (no mean fields are as- 
sumed) with random phases and amplitudes in Fourier 



space. This algorithm ensures conservation of total en- 
ergy and mean fluid density per unit time in the absence 
of charge exchange and external random forcing. Addi- 
tionally, V • B = is satisfied at each time step. Our code 
is massively parallelized using Message Passing Inter- 
face (MPI) libraries to facilitate higher resolution. The 
2D simulations are not only computationally simpler and 
less expensive (compared with the full 3D), they offer 
significantly higher resolutions even on moderate size 
small cluster machines like Beowulf. Since we are inter- 
ested in investgating the inertial range spectra, we restrict 
our simulations to 2D. Our simulations however retain 
all three components of the magnetic field and the back- 
ground field is assumed along the y-direction. This en- 
ables us to treat the magnetic field appropriately The ini- 
tial isotropic turbulent spectrum of fluctuations is chosen 
to be close to kr 2 (where k is characteristic mode) with 
random phases in both x and y directions. The choice of 
such (or even a flatter than -2) spectrum does not influ- 
ence the dynamical evolution as the final state in our sim- 
ulations progresses towards fully developed turbulence. 
While the turbulence code is evolved with time steps 
resolved self-consistently by the coupled fluid motions, 
the nonlinear interaction time scales associated with the 
plasma 1 / (k • U p (k) ) and the neutral 1 / (k • V„ (k) ) fluids 
can obviously be disparate. Accordingly, turbulent trans- 
port of energy in the plasma and the neutral ISM fluids 
takes place on distinctively separate time scales. 

Spectral transfer in partially ionized fluid turbulence 
progresses under the action of nonlinear interactions 
as well as charge exchange sources. Energy cascades 
amongst turbulent eddies of various scale sizes and be- 
tween the plasma and the neutral fluids. In a freely decay- 
ing case, plasma and neutral fluids evolve under the influ- 
ence of charge exchange forces which dramatically affect 
the energy cascades in the inertial range. This is evident 
from Fig. ([TJ where Kolmogorov-like |5|, 14, |6fl fully de- 
veloped turbulent spectra in the inertial range are shown 
respectively for the coupled neutral and plasma fluids. 
The neutral fluid energy spectrum in the inertial regime 
for the coupled system exhibits a k~ a spectrum, where 
the spectral index a « 3.2. On the other hand, the spec- 
tral index for plasma magnetic and kinetic energy spec- 
tra for the coupled system is a « 2.4. The inertial range 
spectral indices for the coupled plasma-neutral system in 
our simulations show a significant deviation from their 
corresponding uncoupled analogues which are respec- 
tively —3 and —1.7 for the neutral and plasma fluids. The 
spectral indices observed in our simulations can be fol- 
lowed from Kolmogorov-like phenomenology JH 0, 0], 
as described in the following. 

The typical nonlinear interaction time-scale (t„/) in or- 
dinary (i.e. uncoupled) plasma and/or neutral turbulence 
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where or V( is the velocity of turbulent eddies. In 
the presence of charge exchange interactions, the ordi- 
nary nonlinear interaction time-scale of fluid turbulence 
is modified by a factor k ce /k such that the new nonlin- 
ear interaction time-scale (Ttvl) in the partially ionized 
heliospheric turbulence is now 



k kv k 
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On using the fact that k ce is typically larger than k (or k Cl 
the characteristic turbulent mode, as defined elsewhere 
in the paper), i.e. k ce /k > 1 in the heliosphere J3>H,[I1], 
the new nonlinear time is k ce /k times bigger than the 
old nonlinear time i.e. Tnl ~ (kee/k)T„i- This enhanced 
nonlinear interaction time in the partially ionized plasma 
is likely to prolong turbulent energy cascade rates. It is 
because of this enhanced or prolonged interaction time 
that a relatively large spectral transfer of turbulent modes 
tends to steepen the inertial range turbulent spectra in 
both plasma and neutral fluids. By extending the above 
phenomenological analysis, one can deduce exact (ana- 
lytic) spectral indices of the inertial range decaying tur- 
bulent spectra, as follows. The new nonlinear interaction 
time-scale of coupled plasma-neutral turbulence can be 
rearranged as 
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where x ce ~ (k ce Vk) 1 represents the charge exchange 
time scale. The energy dissipation rate (e) associated 



with the coupled plasma-neutral system can be deter- 
mined from £ ~ E^/tnl, which leads to 



£ ~ 



Vk 
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where E k is turbulent energy per unit mode. According to 
the Kolmogorov theory, the spectral cascades are local in 
fc-space and the inertial range energy spectrum depends 
upon the energy dissipation rates and the characteristic 
turbulent modes, such that E k ~ e r k^ . Upon substitution 
of the above quantities and equating the power of identi- 
cal bases, one obtains 

£ 2 /V 7 / 3 



(7) 



for the plasma spectrum (the forward cascade inertial 
range). The spectral index associated with this spec- 
trum, i.e. 7/3 w 2.33, is consistent with the plasma spec- 
trum observed in the (coupled plasma-neutral) simula- 
tions (see Fig. (Q}). Similar arguments in the context of 
neutral fluids, when coupled with the plasma fluid in the 
heliosphere, lead to the energy dissipation rates 
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This further yields the forward cascade (neutral) energy 
spectrum 

E k ~e 2 l 3 k- n l* (9) 
which is close to the simulation result shown in Fig. (Q~|). 

4. CHARGE EXCHANGE PROCESSES 

A key issue finally is to understand what role charge 
exchange modes play in coupled plasma-neutral helio- 
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FIGURE 2. Charge exchange spectra 



spheric turbulence as they result essentially from the 
nonlinear charge exchange interactions. To address this 
issue, we plot charge exchange sources associated with 
the momentum and energy equations, Eqs. (Q~|) & (O, in 
Fig. (12). It appears from Fig. (|2} that spectral energy is 
transferred predominantly at the larger scales by means 
of charge exchange mode coupling processes. The latter 
couples the large-scales, or smaller than k ce modes, effi- 
ciently to low-A: turbulent modes. It is primarily because 
of this k ce — k mode-coupling in the smaller k part of the 
spectrum of the coupled plasma-neutral turbulence, that 
energy is pumped efficiently at the lower k inertial range 
turbulent modes. The efficient coupling of the Fourier 
modes at low k's further enhances the nonlinear eddy 
time-scales associated with the coupled plasma-neutral 
turbulence system which is consistent with the scaling 
tNL ~ (k C e/k)Tni, where k ce /k > 1, as described above. 
This consequently leads to the steepening of the inertial 
range spectra observed in Fig. £[). By contrast, higher 
k modes, far from the energy cascade inertial range, are 
notably inefficient in transferring energy and momentum 
via charge exchange mode coupling interactions and are 
damped by small-scale dissipative processes in the cou- 
pled plasma-neutral heliospheric turbulence. 

5. CONCLUSIONS 



change modes modify the heliospheric turbulence cas- 
cades dramatically by enhancing nonlinear interaction 
time-scales on large scales. Thus on scales I > t ce , the 
coupled plasma system evolves differently than the un- 
coupled system where large-scale turbulent fluctuations 
are strongly correlated with charge-exchange modes and 
they efficiently behave as driven (by charge exchange) 
energy containing modes of heliospheric turbulence. By 
contrast, small scale turbulent fluctuations are unaffected 
by charge exchange modes which evolve like the uncou- 
pled system as the latter becomes less important near the 
larger k part of the turbulent spectrum. The neutral fluid, 
under the action of charge exchange, tends to enhance 
the cascade rates by isotropizing the turbulence on a rel- 
atively long time scale. This tends to modify the charac- 
teristics of heliospheric turbulence which can be signifi- 
cantly different from the Kolmogorov phenomenology of 
fully developed turbulence. We believe that, it is because 
of this enhanced nonlinear eddy interaction time, that a 
large spectral transfer of turbulent energy tends to smear 
the current sheets in the magnetic field fluctuations and 
further cascade energy to the lower Fourier modes in the 
inertial range turbulent spectra. Consequently it leads to a 
steeper power spectrum. It is to be noted that the present 
model does not consider an external driving mechanism, 
hence the turbulence is freely decaying. Driven turbu- 
lence, such as due to large scale external forcing, e.g. su- 
pernova explosion, may force turbulence at larger scales. 
This can modify the cascade dynamics in a manner usu- 
ally described by dual cascade process. 
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We have developed a self-consistent fluid model to de- 
scribe nonlinear turbulent processes in a partially ionized 
and magnetized heliospheric gas. The charge exchange 
interactions couple the plasma and the neutral fluids by 
exciting a characteristic charge exchange coupling mode 
k ce , which is different from the characteristic turbulent 
mode k of the coupled system. One of the most impor- 
tant points to emerge from our studies is that charge ex- 



